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It is shown how Luttinger liquids may be studied using sea-bosons. The main advantage 
of the sea-boson method is its ability to provide information about short- wavelength physics 
in addition to the asymptotics and is naturally generalisable to more than one dimension. In 
this article, we solve the Luttinger model and the Calogero-Sutherland model, the latter in 
the weak-coupling limit. The anomalous exponent we obtain in the former case is identical to 
the one obtained by Mattis and Lieb. We also apply this method to solve the two-dimensional 
analog of the Luttinger model and show that the system is a Landau Fermi liquid. Then we 
solve the model of spinless fermions in one-dimension with long-range (gauge) interactions 
and map the Wigner crystal phase of the system. 

I. INTRODUCTION 

In this article, we show how Luttinger liquids may be obtained from sea-bosons. This rectifies some serious 
errors of judgement of previous works^. In addition, we are able to solve for the one-particle properties of 
various systems highlighted in the abstract. Specifically, we are able to write down a closed formula for the 
momentum distribution of the Wigner crystal mentioned in the abstract that does not involve arbitrarily 
chosen momentum cutoffs. Thus our formalism is able to provide information about short-wavelength physics 
as well as the asymptotics. The subject of fermions in one-dimension is vast and we shall only provide some 
representative references. Many of the more important references on this subject have already been included 
in our earlier work^. Some others are by Haldane^, Schulz, Cuniberti and Picri^ and Schulz^. 

II. THE SEA-BOSON METHOD 

In this article, we lay down the prescription for computing the momentum distribution of Fermi systems 
using the sea-boson method. We shall be quite cryptic at times since by now the terminology is well-known^ 
In our earlier work^, we made it quite clear that the sea-displacement annihilation operator has to be 
defined as follows. 

nf(k-q/2)(l-nf(k + q/2)) ^ 
^k(q) = ,„ c;^-q/2Ck+q/2 (1) 

This definition is ambiguous due to the square root of the number operator in the denominator, but leaving 
that aside (see the earlier work'^ for more details), it is clear that we may use this to rewrite the number 
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operator^ as follows (where npO^) = 6{kF — |k|) is the momentum distribution of the free Fermi gas), 

Tik = nir(k) + (1 - nir(k)) (2) 

Here, 

< = i - E 4+q/2 (q)^k+q/2 (q) (3) 
=E4-q/2(q)^k-q/2(q) (4) 

and i = N/N'^. Here N'^ = X^k "^'C') ^'^'^ — ^kCk- From Eq.( 2) it is clear that in general we should 
expect a discontinuity in the momentum distribution at |k| = kp since in general, {n^^) ^ {n^j,}- This 
means that we may regard Eq.( 2) as a proof of the Luttinger theorem. This is true no matter how exactly 
or approximately we have solved our equations. It seems therefore that Landau Fermi liquids are generic. 
But there is a way in which this theorem can be violated. The only way in which Luttinger liquids can 
arise for arbitrarily weak repulsion is if the sums over q in Eq.( 4) and Eq.( 3) diverge in such a way as to 
compensate for the individual terms in the sum being vanishingly small. Indeed, Eq.( 2) suggests a natural 
classification of spinless one component Fermi liquids. We may consider the functions (n^) and (n^) to be 
well defined and nonzero in general, (possibly with infinite slope at |k| = kp) for all k after factoring out 
the discontinuous nF{\s.) and 1 — ni?(k). 
Defn : If we have, 

(i) (^kj?) ^ ('^ky) • L8,ndau Fermi liquid. 

(ii) ("-kr) ~ ("kj!.) • Rcsidueless Fermi system. 

(iii) (^kjj) < (fT-^p) ■ Inverse Landau Fermi liquid. 

It seems that the third possibility has not been investigated at all, and with good reason. Later in this 
article, we argue that there is no room for such a possibility, the most exotic ones are Luttinger liquids (or 
Wigner crystals). 

The operator in Eq.( 1) is not an exact boson but we may treat it as such and impose canonical boson 
commutation rules (this time we include spin for the sake of generality) : [Ak<j(q(T ),Aj^^(q(T )] = nir(k — 
q/2)(l — npCk + q/2)). However now it appears that we must pay a price for this luxury. The reason is 
quite subtle. Treating this object as an exact boson and using the formulas for the momentum distribution 
written down in the main text of our first work leads to logarithimic divergences. In order to tame these 
divergences we have to adopt the following prescription. 

(nk.) = n^(k) {n<J + (1 - n^(k)) (n>,) (5) 

(6) 



exp[2S%{\ia)] 
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exp[2S'i{ka)] 



(7) 



By definition, we have, 

S%(ka) = ^(Gol4+q/2<,(qai)ak+q/2o(qai)|Go)(l - n^(k + q)) (8) 



qcri 



5^(k<7) = ^(Go|at_q/2.i(q'^)«k-q/2.i (qa)|Go)nf (k - q) (9) 

We define, Ai^^{qa') = nir(k — q/2)(l — nF(k + q/2))ako(q(7'). Here |Go) is the ground state of the full 
hamiltonian with the the sea- displacement operators being treated as canonical bosons with the commutation 
rule [ak<7(qcr ), oiko-(*l^ )] = 1 i ^^'^ other commutators involving any two of these operators are zero. Also 
it involves ignoring the presence of the square root making the Fermi bilinear a simple linear combination. 
If q 7^ then, 

4+q/2oCk-q/2o' = Ao(-qa') + (qa) (10) 

+ Yl 4-qi/2ai(qi^)^k-qi/2oi(qi<^') - Yl 4+q,/2o'(^l'^l)^k+qi/2o(qiCTl) (11) 
qio-i qio-i 

From Eqn.( 11) it is clear that we must also have, < n^^ >= 1 — S'B(kc7) and < n^^ >= ^^(kcT). where, 

5s(ka) = ^(G|4^q/2o(q^i)«k+q/2o(q<Ti)|G)(l - n^(k + q)) (12) 

qo-1 

^A(ka) = ^(G|4_^/2^^(qa)ak_q/2oi(q^)|G)n^(k- q) (13) 
qoi 

Here |G) is the ground state of the full hamiltonian obtained by treating the operators Akcr(qo" ) more 
carefully and also by retaining the square root of the number operator and so on. All this is needed 
presumably, since Luttingcr liquids being nonidcal. their momentum distributions fluctuate and one must 
solve the system sclf-consistently. A simple exercise illustrates this point. Take for example the fluctuation 
in the number operator A^(k, k) =< > — < rik >^. From idempotence we may write, A^(k, k) =< 
Uk > (1— < nk >). For a nonitncracting system at zero temperature this quantity is zero. For a Landau 
Fermi liquid with quasiparticlc residue close to unity, this quantity is small. However for a Luttingcr liquid 
this quantity and is large equal to the (square of the) mean itself, since in the vicinity of |k| = kp for a 
Luttingcr liquid, < n^p >« 1/2. This points to the importance of including fluctuations in the momentum 
distribution. In our work this is done implicitly by resumming the infrared divergences. However we have 
found that the ofl'-diagonal counterpart namely, A''(k, k ) =< rik n^' > — < nk >< n^^' > is vanishingly 
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small in the thermodynamic limit. A systematic mathematically rigorous approach is beyond the scope of 
this article, but the claim is the prescriptions in Eq.( 6) and Eq.( 7) capture all these essential features. In 
order to make this plausible, we first note the following feature. Since < S*^ ^ < oo, we have < n^'^ < 1, 
as required for fermions. Also in the case of Landau Fermi liquid for weak coupling, < S*^ ^(kcr) "C 1, this 
leads to the simple formulas presented above with \G) = \Go) and -5^ ^ = Sa,b- In the case of a Luttinger 
liquid, the quantities bO^'^) diverge logarithmically in the vicinity of |k| = kp- The use of Eq.( 6) gives 
us the following momentum distribution of a Luttinger liquid in the case of the Luttinger hamiltonian and 
also in the case of the Calogero-Sutherland hamiltonian( here A is a momentum cutoff, see below ). 



nk = 



e{A-\\k\-kF\ 



+ em - kF\ - A)e{kF - \k\) (14) 

The anomalous exponent 7 is characteristic of the model used. Thus we may now proceed to compute the 
exponent in either case. 

The other question worth answering is the one posed in the beginning, namely is there such a thing as an 
"Inverse Landau Fermi liquid" ? An examination of Eq.( 6) and Eq.( 7) tells us that the answer is "No". 
This is because we find there that (n^^^) > (n^^^). Thus we have accounted for Landau Fermi Hquids, 
Luttinger liquids (and Wigner crystals) and there is no such thing as an Inverse Landau Fermi liquid. 



III. THE LUTTINGER MODEL 
IN THE SEA-BOSON LANGUAGE 

Here we write down the spinless Luttinger model in the sea-displacement language. This is then solved by 
a straightforward diagonalization technique and the momentum distribution is calculated. The hamiltonian 
is as follows. 

^ = E (^) 4mkiQ) + [Siq)Si-q) + S\-q)SHq)] (15) 

kq ^ q 

We have tried to retain a close similarity with the notation in Mattis and Lieb^. Here, the interactions 
couple only the off-diagonal terms, since we anticipate that these are responsible for breaking Fermi liquid 
behavior. Here S{q) = J2k^k{q) C{k,q) and C{k,q) = 6'(A - {kp - \k + q/2\\)e{A - {kp - \k - q/2\\). The 
cutoff A is needed in order to ensure that only electrons in a shell of [ikp — A, ±kF + A] interact and the 
rest are free. We also postulate that v{q) = v{O)0{2A — \q\). In other words, the delta-function repulsion, 
but with \q\ restricted to be small enough. If g > then we find that the sum over k in the definition of S{q) 
is in the region k G [kp — A,kF + A]. Thus this corresponds to right movers or pi in the notation of Lieb 
and Mattis. Similarly, S{—q) involves the sum peaked at k = —kp and corresponds to left movers. Thus 
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this term couples the right-movers to the left-movers and is responsible for breaking Fermi-liquid behaviour. 
In the absence of the 6{2h. — \q\) we find due to the identity below and the repulsion-attraction duality^ 
(more commonly called back scattering) the interaction term vanishes identically ( 'Y^qS{q)S{—q) = ). 
Thus we must really make sure that v{k — k') can be ignored in comparison with v{q). This is possible if 
we include 6{2A - \q\) since then 0{2A - \k- k'\) 0{2A - |2A;f|) = 0. We shall not provide the details of 
the diagonalisation which is by now quite well-known to the readers. The dispersion of the collective mode 
is (jJc{q) = Veff\q\ = vf (^1 — [Xv{0)/ttvf]'^) ^ \q\. Thus we arrive at the following formula for the sea-boson 
number operator. 



< Go\al{q)ak{q)\Go >= 



27r 



\v{0) 



C{k,q) 



L 2Veff{vF + Veff) 



V T^VF J 



- 1 



(16) 



Thus wc sec the emergence of the sorts of exponents found by Lieb and Mattis^. Now we calculate S\{k) 
S%{k)=So{k) as follows. 

So{k) = Y^ < Go\al_^j^{q)ak-q/2{q)\Go > e{kF - \k - q\) 



L 
2^ 



dq 



2ti 



2L 



9iA - \kF - \k\\)e{A - \kF -\k- qWmkp -\k- q\) 



( MO) 

\ TTVf 



= -9{A-\kF-\k\ 



V T^VF ) 



Ln 



A 



\kF - \k\ 



The momentum distribution is given by (using Eq.( 6) Eq.( 7) when — fci?] < A), 



1 



rik 



1 + sgn{kF — \k\) 



\kF — \k\ 



with the exact exponent completely identical to the one in Mattis and Lieb, 

_ 1 

2\ 2 



7 = 



\ ttvf J 



1 



They^ set = 1 which we have displayed explicitly here. 



IV. THE CALOGERO SUTHERLAND MODEL 



(17) 



(18) 



(19) 



Now we would like to see if other systems can also be Luttinger liquids, like the Calogcro-Sutherland model 
(CSM). If one uses the the form Vq = —n(3{(} — l)|g|/(2m) (which is the Fourier transform of the inverse 
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square interaction) in the formulas in our earlier work, we find that the integrals are finite at |fc| = kp and 
this means we have a Landau Fermi liquid for small f3{l3 — 1). This result contradicts the exact solution via 
Jack polynomials^ which shows unequivocally that the CSM is a Luttinger liquid. The reason for this fallacy 
is the repulsion attraction duality explained in our earlier work^. The full hamiltonian consists of a kinetic 
plus exchange part and a correlation part. The correlation part in terms of the sea-bosons that takes into 
account the repulsion attraction duality has already been written down^. We try and solve the model when 
< (3{f3 — 1) <C 1. In general we would have to diagonalize the nonserparable hamiltonian given below. 

^ = E (^) + E E '2!'"'' iM-Q) + 4mAk' w + 4 (-g)i (20) 

feg k^k' 1 

In the weak coupling case we may write, retaining only the terms that break Fermi liquid behaviour. 




Using ideas given earlier, we may read off the anomalous exponent 7 w A^t;^(0)/27r^t;|, = — l))^/8. 
Clearly, for < /3(/3 — 1) <C 1 we have 7 <C 1. But for /? = 3 the anomalous exponent deduced using the 
Mattis Lieb formula is imaginary. This means we have to include forward scattering terms (that is Vq in 
addition to V^_fc' ) that are known to make the exponent real no matter how strong the coupling is. The 
solution of the CSM using Jack polynomials^ tells us that 7 < 1 for /? = 3. ^ This means we really have to 
be more careful when dealing with systems that have |Aw(0)| > ttvf- Unfortunately, we have not succeeded 
in reproducing the results obtained via Jack polynomials even qualitatively, apart from being reasonably 
sure that the system is in fact a Luttinger liquid. This is due to the non-separable nature of the hamiltonian 
(Eq.( 20)) in the sea-boson language brought about by the need to invoke repulsion-attraction duality^. 
We also note that recently Liu^ has shown how to compute the asymptotics of the one-particle correlation 
functions in various dimensions for various systems using a powerful method known as eigenfunctional theory. 

V. THE LUTTINGER MODEL IN TWO DIMENSIONS 

We solve the analog of the Luttinger model in two space dimensions and show that ground state of the 
system is a Landau Fermi liquid. Many authors have in the past addressed this issue. Notable among those 
are works by Anderson^, Baskaran^" and Bares and Wen^^. In the present case, the hamiltonian in the 
sea-boson language is given by, 

^ = E (^) 4(q)^k(q) + E 1^ Y^^M-ci) + Ai{c^)] [A^' (q) + 4 (-q)] (22) 

kq ^ ^ qT^O kk' 



^This requires quite a bit of effort to show even with the propagator written down expHcitly, thanks to D. Sen for 
showing me how. 
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This hamiltonian describes a self-interacting Fermi gas provided we assume that Uq = ^ 9{A — |q|) where 
A <^ kp here < <C 1 is a dimensionless parameter. We may solve for the boson occupation number as 
follows. 

4(q)^.(q))4E7f^..^(-ci) (23) 

where Ak(-q) = nF(k - q/2)(l - nF(k + q/2)). 

„-2( 1 V- nF(k-q/2)-nF(k + q/2) 

k ^ ' ni ' 

r« , A - 1 4- ^1 ^F(k + q/2)-nf(k-q/2) 

k m 

As argued in earlier works, we have to interpret sum over modes with care so as to not lose the particle-hole 

mode, the collective mode being obvious. The is particularly important in two dimensions where we expect 

both to be present. Thus the sum over modes is defined as follows. 

. X du W(q,Lj) f(q, w) , , 

i Jo dujW{ci,u)) 

Here the weight function is given by, 

W (ci, uj) =-Im( — ] (27) 

In our earlier work, we had suggested that in the above formula we have to use a dielectric function that 
is sensitive to significant qualitative changes in one-particle properties. The simple RPA dielectric function 
does not possess qualities we expect from a Wigner crystal. Thus we shall have to derive a new dielectric 
function using the localised basis rather than the plane-wave basis. In the case of the Fermi gas in two 
dimensions we find that the system is a Landau Fermi liquid and there is no need to use a better dielectric 
function, the simple RPA sufiices. The momentum distirbution is always given by, 

■^^g-2S«(k)l 



< nk >= i 



ufO^) + ^ [l - e-^sIKkjJ (1 _ „^(k)) (28) 

S'iik) = J2 (<_q/2(q)^k-q/2(q)) (29) 
q 

S%{k) = ^ (^Uq/2(q)A+q/2(q)) (30) 
q 

The computation of the boson occupation number ^Aj^(q)Ak(q)^ is the key to evaluating one-particle 
properties. 

In the ease of a Fermi gas in two dimensions with short-range repulsive interactions we may use the simple 
RPA dielectric function. The integrals are somewhat complicated since in two dimensions, the angular parts 
are very troublesome unlike in three dimensions. Therefore we copy the results first derived by Stern^^. 
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r / X 1 , mkFVq q 



q muj 
2kF kFq 



q rriLO \ 

— 1 1 

2kF kFq J 



^RPA 



where, 



mkFVq f 

2irq ^ 



C± = sgn 





q 


mui 


2 


2 

-D+ 


1 - 


q mui 
2kF kFq_ 


2" 


1 - 


2kF 


kFq. 









} (31) 



(32) 



q 




, D± = 0, 


q 






_l_ muj 






_2kF 


kFq. 


2kF 


kFq 



> 1 



9 ^(-q,w) 



mkF ^mC- 



2'Kq kFq 



C± = 0, D± = 1, 



q mwV 
2kF kFq J 



q _|_ mu 
2kF kFq 



< 1 



(33) 



(34) 



\2kF kFq 



kFq 



\2kF kFq) 



\2kF kFq)' 
(35) 



4(q)A.(q)) = ^^^'^ '^^^''^^ (ffS)^ ^^^-^'"^ 



W{q,uj) = 



iq,i^) 



^RPAiq^^) + ^kpAiq^^) 

POO 

Z{q) = doj W{q,uj) 
Jo 



In general we have, 



1 f°° 1 

SaO^) ^ jTTyl I '^l 1 / '^'^ W{q,uj) fA{k,q,w) g'^{-ci,uj) 

[Iny Jo Z[q) Jo 

1 f°° I f°° 

^"^^^'^ " (2^ Jo ^ ^{q) Jo ^(9'^^) fB{k,q,U!)g^{-q,uj) 



(36) 



(37) 



(38) 



(39) 



The rest of the details are relegated to Appendix A. The collective mode occurs when Im[e] = 0, that is, for 
small enough q. This means that we have to treat this separately. 

1 

WTJo "'^^ Z{q) 



1 f°° 

dq q dui W{q,uj) fA{k,q,u;) 3^(-q,a;) 



+ 



1 f°° 

j^;^ J dq q /^(A;, g, Wc) 5^(-q, Wc) 



(40) 



dq q 



(27r)2 Jo ^ ^ Z{q) Jo 



duj W{q,uj) fBik,q,uj)g^{-q,u}) 
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1 f°° 

+J^;^ J dq q fB{k,q,uJc)g'^{-ci,LJc) (41) 

Here it is implicit that in W we assume that Im[e] ^ 0. The dispersion of the collective mode may be found 
using Mathematical'^ . It is given below. 

{Tr^q'^+rmTq'^Vg + klm''v^)^ 

This dispersion is real and positive for all q and for all Vq > 0. Thus in the small q limit, where using just 
the RPA dielectric function is justified and is also the limit where the close to the Fermi surface features of 
the momentum distribution is given exactly, we are justified in retaining only the coherent part. Thus we 
may write, 



1 f°° 

^AiK^j^:^ J dqq fA{k,q,Wc) g^{-ii,uJc 



(43) 



1 

^^(k) « ^ dq q fBik,q,iOe)9H-n,^c) (44) 

To determine whether or not Fermi liquid theory breaks down, we have to compute, 

1 f°° 

^A{kF)^j^^j dq q fA{kF,q,u)c) 9'^{-<l,0Jc) (45) 
1 f°° 

S%{kF) ^ jTT-r^ dq q fB{kF,q,uJc)g'^i-ci,Wc} (46) 

l^^J Jo 

If S^ikp), SB{kF) < oo then the ground state is a Landau Fermi liquid. If S^ikp) = >S'g(fci?) = oo then the 
system is a non- Fermi liquid. For small q if we set = v^f f q we have, 

fAikp, q, Wc) ~ fnikp, q, Wc) ~ 1/?^ (47) 

Also, 

g^{-(l,u;c) ^ q (48) 

Thus the integrals in Eq.( 45) and Eq.( 46) are infrared finite. This means that S%{kF), S%{kF) < oo and 
the system is a Landau Fermi liquid. The details of the momentum distribution can be worked out but are 
not very important. A closed formula for the quasiparticle residue may be written down as shown below. 

Zf = ^ (49) 
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VI. ONE DIMENSIONAL SYSTEM WITH LONG-RANGE INTERACTIONS 



In this section, we consider electrons on a circle interacting via a two-body attractive long-range interaction 
with strength proportional to the separation between the electrons. In this case, we expect the system to be 
a Wigner crystal since the interaction is long-range and actually increases (in magnitude) with separation 
rather than decreases. This means that the electrons prefer to be as far apart from each other as possible to 
lower the energy leading to a crystalline ground state. Thus we have to be careful about the choice of the 
dielectric function. First, we postulate that Vg = 2e^/ {q a)^ which corresponds to the gauge potential. Here 
a has dimensions of length and > is dimensionless. From the form of this potential, one hopes that we 
need not concern ourselves with the issues that were relevant in the case of the Calogero-Sutherland model 
namely the repulsion attraction duality. Thus we may write as before, 



4(q)^k(q) 



rcLu W{q,co) ^''^^'^^ 9\-ci,co) 
V Z{q) Jq ^ ' c,.,-u!s:a^2 y V "1' ^ 



(w + 



W{q,u:) = 



e^^{q,Lo)+ei^{q,Lo) 



(50) 



/■OO 

Z{q) = doj W{q,uj) 
Jo 



(51) 



1 d 

fif"2(-q,w) = — — e(q,w) 

Vq OLO 



(52) 



As mentioned before, we have to be extra careful in making sure that we choose the right dielectric function. 
The RPA-dielectric function is not likely to suffice since its static structure factor (SSF) does not exhibit 
the features we expect from a Wigner crystal. In particular, we expect 5(2^^) = oo as we shall see soon. 
To convince ourselves of this we ascertain the properties of the RPA dielectric function with long-range 
interactions. 

2- 



m 



{kF + ql2f~(^^y 



{kp-q/2f-[^) _ 
The zero of the above dielectric function gives us the dispersion of the collective modes. 



(53) 



Wc(g) = 



m 



{kF + q/2y-ikF-q/2rexp{-^) 



(54) 



For 1^1 -C kp and Vq = 2e /{aq) we find. 



^c{q)^ -\ - + W-o +0{q ) 



m 



(55) 
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This plasmon-like gap loq = ^\ ^^^f ™ \ in the collective mode is present due to the characteristic 
nature of the potential. But this is also present in the three dimensional electron gas and is not a sign 
of an insulator since the latter is not at high densities. A gap in the one-particle Green function could 
be taken as a sign of insulating bchaviour^^. However, in our approach we are unable to compute the full 
Green function as yet. Thus wc must resort to a more indirect approach. For a Wigner crystal, the SSF 
must exhibit certain singularities. Thus we have to use the generalised- RPA that is sensitive to qualitative 
changes in single-particle properties. The new dielectric function will involve the full momentum distribution 
which has to be determined sclf-consistcntly using the above sea-boson equations. In our earlier work we 
suggested that the new dielectric function should also involve fluctutations in the momentum distribution, 
however it now appears that that is fortunately not needed. The number-number correlation function is 
vanishingly small in the thermodynamic limit as shown in another preprint and this means we may simply 
write, 

e(q,c.) = l + ^E /T"!7^^ (56) 

and the momentum distribution is determined self-consistently using the sea-boson equations (Eq.( 28)). 
This is too difficult to solve analytically and hence we have to resort to a numerical solution. In order to 
simplify proceedings even further, we use only the collective mode. The particle-hole mode which is due to 
a nonzero Im[e] is needed if one is interested in features of the momentum distribution away from the Fermi 
surface more accurately. However we shall hope that this is given not too badly even at these regions far 
from the Fermi points. 



4(q)Ak(q)\ = - ^''^ 5'(-q,Wc) 



9 ^(-q,Wc) 



Trq 



.^c - {'"fQ + ujI - [vpq - Cqf 



(58) 



(59) 



^q ('^c(g) + ^ - eg)2 

Sl{k) = i V (^'""f + '^)^ ,2(-q,a;e) (61) 

To proceed further, we have to ascertain the nature of the collective modes ujc- If we use the RPA-dielectric 
function, we find a constant dispersion (plasmon) for small \q\. However we have found that this choice is 
inconsistent since if we use the momentum distribution obtained from this to solve for the dielectric function 
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and recompute the collective mode we obtain a completely different answer namely : uJcio) = Vs\q\- Therefore 
it is critical that we get the dispersion right. It appears then that we have to use the form given in the 
appendix which is not easy to simplify. A systematic approach for obtaining the dispersion of the collective 
modes has been suggested by Sen and Baskaran^°. Since the plane- wave basis is not appropriate for deriving 
a formula for the dielectric function of a Wigner crystal, we shall follow this approach. First, we would 
like ascertain the lattice structure in the small a limit. In this limit, the potential energy dominates over 
the kinetic energy. If we assume that the electrons are all on a circle of perimeter L then to minimise the 
potential energy, we have to maximize the separation. This leads to an equally spaced set of lattice points 
with lattice constant Ic such that N l^. = L. Thus we have Ic = 1/po = Tr/kp- Thus we assume that the 
electrons all lie on a circle with equal spacing between them. Therefore we expect the structure factor to 
diverge for a momentum qo = 27r//c = 2kp. From the Bijl-Feynman formula S{q) = eq/Ludq) we may suspect 
that a choice of lUc that vanishes at q = 2kp is needed. The form of the dispersion is given in the appendix. 
For -k/N ^ \q lc\ ^ 27r it seems that ujq w t^o- For \q lc\ ^ tt/TV we have to be more careful. And of 



course we must have < 



for q Ic = ±27r. But since in the thermodynamic limit tt/N « wc may choose 



(hopefully) LOq ^ luq. In Fig. 1 and 2 wc sec the momentum distribution obtained from these formulas has 
been plotted. In fact, we may write down a closed formula for the momentum distribution. 



rife 



-{1 + Exp 



m uio 



k'p 



fc2 



)nF{k) + ^{l-Exp 



m LOq 
^2 _ p 



){l-np{k)) 



(62) 



The striking feature of this momentum distribution is that it is perfectly flat at = kp- In other words, 
not only is the slope zero but all the derivatives of the momentum distribution vanish at |fc| = kp. This is a 
striking prediction. This may be contrasted with the smooth Gaussian function of Gori-Giorgi and Ziesche^^ ( 
Eq.(Bl) in their Appendix B ). But they consider three dimensional systems which may be different from the 
one studied here. One particle spectral functions are accessible to tunneling experiments or angle-resolved 
photoemission spectroscopy (ARPES). A more difficult problem may be to experimentally realise a Id electron 
system with long-range gauge interactions. 

0.7 




1.01 
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FIG. 1. Momentum Distribution of a Wigner Crystal (with zoom) 



^ 0.5 




FIG. 2. Momentum Distribution of a Wigner Crystal (no zoom) 

The formula below for the static structure factor is derived in the appendix. 

1 W(^) _^ 

This may be further simplified in the thermodynamic limit as follows. Consider, 

sin{N x) 



S{x) 



TT X 



Then we may write, 



S{q) = Siq)— + d{q - 2fc^)-^ e"^ 



(63) 



(64) 



(65) 



As we can see here the strucutre factor diverges at \q\ =2kF which means the system is a Wigner crystal ( 
a true Wigner crystal, since the divergence is from a delta- function ). 

VII. CONCLUSIONS 



To conclude, we have finally shown how to reproduce Luttinger liquids using sea-bosons. We have also 
shown how systems with both long-range and short-range interactions may be studied with almost equal ease 
and that too in any number of dimensions. Using the same formalism we may also probe short- wavelength 
physics as illustrated by the momentum distribution of the Wigner crystal that contains no momentum 
cutoffs. Finally we showed that the two dimensional analog of the Luttinger model is a Landau Fermi liquid. 
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IX. APPENDIX A 

F_ 

{u, + ^cos{9)-e,r 



fA{k,q,u}) = de — — — (66) 



fB{k,q,iv)= de ^—k-q ,Q. ^ ^2 



Define, 



u{x\ A,B) = I dx 



{A + B Cos[x]f 



Prom Mathematical'^ we find, 



ArcTanh 



u{x;A,B) = -2A- 



{A-B)Tan[%] 



(A2 - 52^52 _ ^2 



B Sin[x\ 



{B"^ - A^){A + B Cos[x]) 



Integrating by parts we find, 



fA{k,q,uj) 



/■27r 

Jo 



ddsinO 5{^^ " ^ 



2kq +cos{9))u{e;u;-e„^) 



r2TT 

w)= / dO sinO 5{ 
Jo 

This may be rewritten as, 

fA{k,q,ui) = 

Jo 



2kq m 



_ i-2 _ 2 
^ 2kq 



+ cos(6)) u(9; u - €g,—) 
m 



(67) 



(68) 



(69) 

(70) 
(71) 



-L 



de sine s{^^ — ^ 

^ 2kq 



— cos(e)) u(e + tt; w — Cg, — ) 

m 



(72) 



fB{k,q,w) = ( dO sinO — + cos{e)) u{G;io + Cg,—) 

Jo 



2kq 



m 
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de sine 5{ 
2fcg 



^ ~l~ Q. kjp / r\\\ //I \ 

— cos(u)) u(u + Tr;uj + eg, — ) 

m 



(73) 



60 = arccos 



2kq 



(74) 



2kq 



(75) 



fA{k,q,w) = u{eo;uj- eg, ■^) 



0{ 



1.2 _ 1,2 _ „2 1.2 _ 1.2 _ „2 

H-F 1 _|. 1] _ Q( F 'i _ 



2kq 



2kq 



-w(C^o + w - eg, — ) 
m 



fc^J 9_ 

. ^ 2kq 



(76) 



fB{k,q,u}) = u{e'o;u + eg, — ) 



-I- o2 _ i.2 1.2 I „2 _ jL2 

^( + 1) - g( 1„ ^ - 1) 



-u(6'o + 7r;a; + Cg, — ) 
m 



0(1 



k^ -\- fc|^^ 6^( 1 fe'^ + 9^ fef ^ 



(77) 



X. APPENDIX B 



Here we use the approach suggested by Sen^° to derive a formula for the collective modes. The formula 
Eq.( 56) although probably right is not very illuminating, for it is hard to see how the structure factor derived 
from this formula possesses the features we expect namely a divergence at |g| =2kF- Thus we would like to 
derive a formula for the dielectric function where this feature is manifest. To do this we adopt the localised 
basis rather than the plane-wave basis. In real space, the hamiltonian we are studying is written as follows. 

JV-l 

ti\ (78) 



/7= V ^--Vlar - 



2m a? 

i=0 i>J 



Wc assume that particles are on a circle and |x| is the chord length. We would like to compute the dielectric 
function using this model. The density operator in momentum space is, 

N-l 



(79) 



m=0 



Here Ic + Xm is measured along the circumference of the circle (see Fig. 3 below). 
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FIG. 3. Schematic Diagram of Electrons on a Circular Lattice 



The average density is given by, 



Af-l 



m— 



It can be shown that (see below) < >= l/mcjo, independent of the index n. Thus we have, 



1 _ e«9^ 
1 - e 



In those instances where < Pg >7^ the static structure factor is given by, 



TV TV 



1 - e'9 



1 sm2(^) 
e ""^0 = — 4 — e """0 



(80) 



(81) 



(82) 



The rest of the details are as follows. We write Xm{t) « R 9m{t). In terms of the small angles 9i the 
hamiltonian in Eq.( 78) may be written as follows. 



Af-l 



H = 



1 V 



(cos(27r— + e^} - cos{27T— + 9^')f + (sm(27r— + 9^) - sin{27r— 



(83) 



We may expand the above hamiltonian in powers of the angle and retain only the leading terms to arrive at 
the following hamiltonian in the harmonic approximation. 



N-l 



p;2 



where, 



A{n, n ) 



ire 
ALo? 



B{n,n ) = - ^ sgn{sin{-K ^ 



[n — n ) (n — n ) 



)) cos{tt 



(84) 

(85) 
(86) 
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Despite appearences to the contrary, the extremum of the potential is at i„ = 0. Since A> 0, this extremum 
is also a minimum. One has to now compute the various correlation function of the system. The primary 
one of interest is, 

Gii{nt;n't) =< T i„(t) x^>{t') > (87) 

The other is, 

C?2i {nt; nt)=<Tpn (t) {t ) > (88) 

Thus we have, 

i—Gxi{nt;nt) = ^G2i{nt;n't) (89) 

i — G2i(nt; n't) = 5^^^>5{t - i') - 4i ^ A(n, j)(Gii(ni; nt) - Gn{jt; n't')) (90) 
This may be solved by a Fourier transform. 

G,,{nt;nt) = ^Yl^'^^'~'' ^ ^12^'' '=("-"')G,,(g, z,) (91) 



-i(3^ N- 

P Q 



Thus we have. 



Thus, 



i Zp Gil {q, Zp) = —G21 {q, Zp) (92) 

i Zp G2i{q,Zp) = 1 + 4i {A{q) - 1(0)) 611(9, ^p) (93) 

A{q) = Y^A{j)e^'^'^^ (94) 
j 

Gil {q, Zp) =(imzl-^i {A{q) - i(0))) (95) 

Gii{nt;nt) = ^E^'' ^ ^^..(t-t') (^^ zl - A {A{q) - A^y^ (96) 

p 



2i / 2i V 1 - e'(* 



1 + e*« 1 + 



(97) 



^(0) = ? f - « ^ (98) 



Here cq = 7re /4La . The dispersion of the collective mode is then given by, 
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u;, = ^^i^A{0)-A{q)y (99) 

U q = ov q = 2kF then Wq = 0. Notice that unless lc\ <C 'k/N, we have A{q) w for \q\ ^ since 
it is vanishingly small in the thermodynamic limit. Thus we have to take the large A'' limit first in which 
case A{q) w for -k/N <^ \q lc\ <C 27r. Thus we may set Ug w wq with impunity in this case. For a 
more thorough analysis, one has to compute the full dielectric function from the dynamical density-density 
correlation function and use it to compute the full momentum distribution that is accurate even away from 
the Fermi surface. However we shall be content at features close to the Fermi surface. For the static structure 
factor we have to compute the equal time version of the correlation function. We find that Wg is in general, 
complex. This means that the eigenmodes also have a finite lifetime. Thus we have, 

1 /I pOD 



Q 



1 . . piq lain-n ) 



2m N ^ uJc{q) 

From Eq.( 99) it is clear that for n/N 4C \q lc\ < 27r we have Wg w Wq since A{q) w for g' in this region. 
Since in the thermodynamic limit, this is all of q, we shall write, 

1 gil lc{n-n) ^ 
< Xn{t) X^' {t) >= -: > = (5„ „' (101) 
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